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ABSTRACT 


The radar response fo soil moisfure contenf was experimenl-ally defermined 

for each of fhree bare fields wifh considerably different surface roughnesses at 

eight frequencies in the 2-8 GHz band for HH and W polarizations. Analysis of 

the data indicates that the effect of roughness on the radar backscattering coefficient 

can be minimized by proper choice of the radar parameters. If, in addition, sensitivity 

to soil moisture variations and system design constraints are considered, the following 

radar parameters for an operational soil moisture mapper are recommended; frequency = 

4 GHz , angle of Incidence range = 7°-15° and either HH or VV polarization. The 

3 

corresponding sensitivity is about 0.25 dB/.Ol g/cm . 
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1.0 INTRODUCTION 

Soil moist'ure cont'ent' is on import'ani' poromefer in crop yield prediction, 
runoff prediction ond other applications in hydrology, ogricul ture and meteoiology . 
Remote sensing offers o potential means for the determination of the spatial distri" 
bution of moisture content over large areas, within a short time and at a reosonable 
cost. Among the various electromagnetic sensors, radar is the least affected by 
atmospheric conditions. It is also time"of“day independent and has the capability 
of mapping terrain surfaces from satellite altitudes with a spatial resolution compatible 
with the requirements of the obove applications. 

The backscattering coefficient, a°, of an area extended target, such as a 
soil medium, is a function of the soil surface roughness and dielectric properties. 

At microwave frequencies, the dielectric constant of soils is strongly dependent upon 
the soil moisture content [1,2]. Radar backscatter measurements of bore soil surfaces 
indicate a high degree of sensitivity to moisture variations, particularly at angles of 
incidence close to nadir [3,4]. Similar observations have also been noted for 
vegetation covered fields [5,6], although the presence of the vegetation tends to 
reduce the sensitivity of cr° fo soil moisture variations. It was further observed that 
differences in vegetation cover and soil surfoce roughness affected both, the absolute 
value of cr° (for a given moisture content) as well as its sensitivity to moisture vari~ 
ations. For vegetated fields, the correlation coefficient between and soil moisture 
content was highest at an angle of incidence of 10° (from nadir) [5] . 

The present study is a detailed analysis of the effect of roughness on of 
bare fields. The objective is to experimentally determine the "optimum" combination(s) 
of radar parameters (frequency, polarization, angle of incidence range) for mapping 
soil moisture content. An optimum combination of sensor parameters is defined here 
such that cr° of the ground is almost independent of surface roughness while retaining 
an acceptable sensitivity to soil moisture variations. The present study is limited to 
bare ground; the combined effects of surface roughness and vegetation cover will be 
the subject of a future investigation. 


2.0 THE EXPERIMENT 


Radar backscaffer dafa were acquired from ^hree bare fields having appreciably 
differenf surface roughnesses using fhe University of Kansas truck-mounted 2-8 GHz 
Microwave Active Spectrometer (MAS) system [7]. The operational characteristics 
of the rodar system are given in Table 1 . 

N A schematic of the test site is shown in Figure 1. For reference purposes, the 
three fields are designated by S, M and R for smooth, medium rough, and rough, 
respectively. Photographs of the three fields are shown in Figure 2. By inserting a 
4' X 2' X 0.5" metal plate into the soil ond photographing it, the surface profiles 
shown in Figure 3 were obtained. After digitizing each of the profiles at the rate of 
10 points per cm (in the horizontal direction), the rms heights were calculated with 
respect to the mean slope of the ground and found to be 0,88 cm, 2.6 cm, and 4,3 cm 
for fields S, M and R, respectively. In terms of peak-to-peak variations, typical 
values for fields S, M and R were 2 cm, 7.5 cm and 15 cm. The combination of rms 
height and peak-to-peak variation clearly illustrates the wide range of surface rough- 
ness represented by the three fields. The surfaces of fields R and M ore typical of 
plowed and plowed-then-disked surface configurations, but the surface of field S is o 
rather rare case in terms of the majority of terrain surfaces. 

After acquiring one radar data set per field with the fields dry, the fields were 
sprinkler irrigated for 12 hours. Radar data sets were then acquired from the three fields 
on a rotating basis for a period of 13 days. Each radar data set consisted of measure- 
ments of the backscattered return at five angles of incidence (0° (nadir) through 40° in 
10° steos) for HH and VV polarizations (see Table 1) at each of eight frequencies in 
the 2-8 GHz bond. To reduce signol foding, o combination of frequency and spatial 
averaging was used. Frequency averaging was provided by the 450 MHz bandwidth 
of the MAS system and spatial averaging was achieved by measuring the return from 
several different spots on the field. Since the effects of signol fading are most severe 
at nadir [8], more spot 'ally independent meosurements were ocquired at nadir than at 
the higher angles of incidence. The number of measurements at each angle and the 
cciculated 50% and 90% confidence intervals associated with the date at that ongle 
are given in Table 2. 

Simultaneous with each radar dato se% soil somples were collected at each 
of the eight sompling locations 'Figure 1 » of the field under observation. At eoch 
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Table 1: MAS 2“8 System Specifications 


Type: 

Modulating Waveform: 
Center Frequencies: 


FM-CW 

Triangular 

2.75 , 3.25 , 4.75 , 5.25 , 5.75 , 6.25, 
6.75, 7.25 GHz 


FM Sweep: AF 
Transmitter Powen 
IF Frequency: Fjp 
IF Bandwidth: AF|p 
Antennas: 

Height above ground: 
Transmitting antenna diameter: 
Receiving antenna diameter: 
Feeds: 

Effective Two“way beamwidth: 
Incidence angle range: 

Polarization: 


450 MHz 
40 mW 
50 KHz 
6 KHz 


20 m 
91 .5 cm 
91 .5 cm 
Log periodic 

5 . 4° at 2.75 GHz/ 2 . 2® at 7.25 GHz 
0® (n«idir)-80® 

Horizontal transmit-Horizontal receive (HH) 
Vertical transmit-Vertical receive (W) 


Calibration: 

Internal 

External 


Delay line 
Luneberg lens 
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Table 2: Calcula^ed 50% and 90% confidence limits around the 

scattering coefficient along with the number of spacially 
discrete measurements acquired at each angle. 


Angle of 
Incidence 

Number of Specially 
Discrete Measurements 
Collected 

Confidence Limits* (dB) 

50% 

90% 

0 

15 

+0.6 

+1.6 



-0.65 

-2.0 

10 

13 

+0.65 

+1.65 



-0.7 

-2.05 

20 

11 

+0.7 

+1.4 



-1.0 

-1.6 

30 

9 

+0.5 

+1.1 



-0.6 

-1.2 

40 

7 

+0.5 

+1.05 



mm 

-1.2 


* Calculated at 5 GHz. 
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X Indicates Soil Sampling Locations 


Figure 1. Schematic of the test site located at Texas A & M 
University agricultural experiment station. 
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Photographs of bore fields. Soil type is Miller Cloy (49^/. cloy, 35% 
silt, and 16% sand.) (a)Smooth surface/ PMS height = 0.88 cm. 

(b) Medium rough surface/ RMS height = 2.6 cm. (c) Rough surfoce / 
RMS height = 4.3 cm. 


F a • ' 



r75 :> 

mm 
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■ ' 






(b) 


(c) 


Figure 3. 



Photographs showing the method used to obtain the surface roughness 
profile of (a) the smooth field, (b) the medium rough field, and (c) 
the rough field. The scale on the metal plate is two inches per 

division . 
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locafion samples were collecfed from five deplfis: (^1 cm, 1-2 cm, 2-5 cm, 5~9 cm, 
and 9-15 cm. These samples were later processed in the laboratory to determine their 
moisture contents by weight and their bulk densities. Complete listing of the ground- 
truth and radar data were documented In a separate report [9] . 


3.0 DEFINITION OF SOIL MOISTURE CONTENT 


In the case of o homogeneous soil medium (in terms of its volumetric water 
content), there is only one obvious quuntitative definition of moisture content. The 
analysis of data to determine the dependence on moisture content Is straightforward. 
Under natural conditions, however, the soil moisture profile with depth is a constant 
only In the extreme coses of very dry or fully saturated soil. For the in-between coses, 
the shape of the soil moisture profile of a given field is a complex function of a variety 
of soil and environmental parameters including texture profile, bulk density profile, 
temperature profile, rain and irrigation history, and others. Since the dielectric con- 
stant of soil is strongly influenced by its moisture content (Figures 4-6), associated 
with a given moisture profile are two profiles representing the real and imaginary parts 
of the relative dielectric constant, and K 2 respectively. At a given microwave 
frequency, the and <2 profiles specify the attenuation profile in the soil medium 
from which the skin depth can be determined. As an illustration. Figure 7a shows the 
skin depth as a function of moisture content calculoted for a homogeneous medium. 

It is clear from Figure 7a that at a frequency of 10 GHz, for example, the moisture 
content in the top 1 cm can be used os an adequate definition of the moisture content 
responsible for the observed radar return (except for very dry soils). 

At a frequency around 1 .3 GHz, on the other hand, the backscattered power 
from a soil medium with o variable moisture profile would most likely include contri- 
butions from the ali-soil interface os well os from subsurface layers. Thus, in order 
to analyze the scatiering coefficient data in terms of moisture content, surfoce rough- 
ness, soil type, or any other variables of interest, it is necessary to develop a method 
by which the moisture profile can be represented by an "equivalent" moisture content. 
The remainder of this section is devoted to a discussion of several approaches for 
"defining" moisture content. ’ 


8 



FREQUENCY: 1. 1 GH/ 



Ftqure 4. Representative dielectric constant vaiues as a function 
of volumetric water content for sand, loam, and clay 
at 1.3 GHz. 



Figure 5. Representative dielectric constant values as a function 
of volumetric water content for sand. loam, and clay 
at 4. 0 GHz. 


fTREOUENCY: 10.0GhT| 



'In 




Figure 6 Representative dielectric constant values as a function 
of volumetric water content lor sand loam, and clay 
at 10. 0 GHz. 
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Figure 7. (a) Skin dep^h and (b) power reflecHon coefficienf as a funcfion of volumefric 
water content, frequency and soil type. 
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3.1 Moisture Wifhin Fixed Dep^hs 

Soil moisture content in the layer between depths o and b below the surface can 
be calculated either on o dry weight basis (gravimetric): 


mw(a/b) = 


t 

/ 


W (z)- W.(z)dz 
— ' dimensionless 


( 1 ) 


or on a volumetric basis: 




b (W,^.(z)- W^(z)]pjz) dz 


W^(z) 


, g/cm'’ (2) 


where 


z = depth below the surface, cm 

o,b= depths of upper and lower ends of the layer, respectively, cm 
W^(z)= profile of the wet weight of the soil, g 


W ,(z) = profile of the dry weight of the soil, g 

® 3 

p(z) = bulk density profile of the soil, g/cm 

s 3 

p^ = density of water = 1 g/cm 


Since the surface layer is the most influential in terms of the sensor's response, m^ and 
m^ are usually evaluated for a - 0. This method has the advantage of simplicity and 
is probably adequate at frequencies around 10 GHz or higher if b is chosen to be 
around 1cm. Between m^ and m^, the latter is preferred because 1) it incorporates 
the variation in bulk density with depth and hence is independent of soil compaction 
and 2) it is physically a better descriptor of the dielectric constant of the soil since 
it is the total number of water molecules per unit volume of the soil“water mixture 
that determines the dielectric constant [ i). Hence, oil moisture ccnts!if« discussed 
in this report are calculated on a volumetric basis. 
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3.2 Moisfure Wifhin ^he Skin Depth 


In an affempf fo construct the soil moisture profile from multifrequency 
microwave radiometric measurements, Poe [10] defined the moisture content of the 
soil in terms of the water within the skin depth. The method was later adapted by 
Ulaby et al . [ 4 ] to the analysis of radar backscatter data from bare soil . 

After dividing the soil medium into horizontal layers each having a thick- 
ness d, the measured profiles of W (z), W ,(z) and P (z), are used to calculate the 

fh w Q s 

volumetric moisture of the i layer rn^(i), for each i between the surface and the 
desired depth. The top layer is i = 1 . Using the relationships between m^ and 

and (Figure 4-6), K,(i) and K«(i) can be determined from which the attenuation 

. . .th ^ 

coefficient of the I layer can be calculated: 



where A is the wavelength. The skin depth 6 is defined by: 

6 

^(z) dz = 1 neper (4) 

o 


For the discrete cose 6 can be determined by summing the attenuation of the layers from 
the surface downward until the total attenuation Is approximately equal to 1 neper, 

6 is then the product of the number of byers used in the summation by the width d. 
Using 6 as a parameter characterizing the moisture profile of the soil medium under 
investigation, on equivalent moisture content of a homogeneous mediunvm^ ,can be 
defined such that both media have the same skin depth 6 . This equivalent moisture 
content can be obtained from 6 using Figure 7a. 
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Although the skin depth method takes the attenuation of the wove into 
account, it does not consider the reflection at the air-soil interface and at the 
interfaces between the various layers. 

3.3 Moisture Content Estimated by o Coherent Model 

An equivalent moisture content of an inhomogeneous soil medium (with depth) 
can be defined os the moisture content of o homogeneous soil medium whose power 
reflection coefficient at the air-scil interface is equal to the power reflection coef- 
ficient of the inhomogeneous medium. The power reflection coefficient of o homogeneous 
medium at nadir is shown in Figure 7b os a function of moisture content, m^. 

The reflection coefficient of a soil medium with varying dielectric constant can 
be calculated from the solution of Maxwell's equations in o plane-stratified dielectric 
medium. Casey [11] formulated the solution to this problem in terms of Hill's functions 
which con be evaluated numerically for a given dielectric profile. Since the method 
treats the cose of a coherent signal incident upon o dielectric medium, it is referred 
to herein os the Coherent Model, as distinguished from the Incoherent Model discussed 
in the next section. 

The Coherent Model incorporates both attenuation and phase shift in the soil 
medium which causes the calculated reflection coefficient to be very sensitive to the 
derivative of the dielectric constant with depth. Hence, unless the soil moisture 
variation with depth is small over on interval of one wavelength, correct evaluation 
of the reflection coefficient necessitates accurate knowledge of the moisture profile 
at depth intervals of the order of 0.1 \ where \ is the wavelength in the medium. 
As on example, at a frequency of 2 GHz, in a medium with a relative dielectric 
constant = 25 is 3 cm. From o practical standpoint, collecting soil samples smaller 
than 1 cm in depth is difficult and time consuming. 

3.4 Moisture Content Estimoted by on Incoherent Model 

Although the wove transmitted by the radar is approximately coherent in nature, 
the wove bockscottered by o target composed of a large number of random scatterers, 
such as a soil medium, is usually described statistically in o manner similar to the 
statistics of random noise [12]. This is a result of spatial or frequency averaging (or 
both). This noise”! ike (incoherent) description of the bockscottered power suggests 
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fhaf the moisfure confenf represenfafion of a soil medium should be bosc-d on fhe 
power reflecHon coefficient calculated for an incoherent signal. For a plane 
surface, the incoherent power reflection coefficient R can be readily calculated 
from the emissivity of the medium f os follows; 

R = 1 - f (5) 

If the medium is homogeneous, the coherent and incoherent power reflection coefficients 
are the same. Burke and Paris [13] constructed a numerical procedure for calculating 
the brightness tempeioture of a layered medium in terms of the physical temperature 
and dielectric constant profiles of the soil medium. If the physicol temperature profile 
is assumed constant, their model reduces to an evaluation of € . 

The fol'owing formulation is adapted from Burke and Paris [13]. With the 
soil medium divided into n homogeneous layers, the emissivity e at nadir is given by: 


e 


(1 - RJ 



z 


i=2 


r. p. n Q. + II 
' ' ■ j=l J / j=1 


Q. 

J 


( 6 ) 


where 


R^ = power reflection coefficient at air-soil interface 

I = layer index (i=l is surface layer and i=n is a semi-infinite 
homogeneous layer) 

r. = 1 - i/L. 

I I 

P. = 1 +R./L. 

I I I 

Q. = (1-R.V'L. 


L 

R 

a 

d 


= exp (a. d.) 

= power reflection coefficient at the interface between the i and i-*-l layers. 
= attenuation coefficient of layer i, nepers/cm 
= width of i^^ layer, cm. 
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From a physical standpoinf', the Incoherent Model is limited to coses where tho 
emission is primarily a consequence of absorption in which cose volume scattering can 
be neglected. This condition is satisfied provided the moisture content in the soil is 
not very small . 

3,5 Comparison of Various Soil Moisture Estimates 

Figure 8 compares five soil moisture estimates at three frequencies for the 

14 profiles obtained from the medium rough field. The five moisture estimates are; 

m^ = moisture content in top 1 cm 

m 2 = moisture content in top 2 cm 

m^ = skin depth equivalent moisture content 

m = Coherent Model equivalent moisture content 
c 

mj = Incoherent Model equivalent moisture content 
The profiles as a function of depth are shown in Figure 9. At the lowest frequency 
(2.75 GHz), for profiles 2, 3, and 7, whose top layers are very moist, all five 
estimates are quite similar. For the driest condition (profile 1) except for the effective 
moisture in the skin depth, the remaining four estimates are similar. As the profiles 
start drying out at the surface the estimates start differing quite appreciably. These 
patterns are similar at 5.25 GHz and 7,25 GHz. 

At the lowest frequency for which the skin depth is comparatively large, the 
skin depth estimator yields the highest estimate. As the frequency is increased the 
skin depth decreases resulting in a m^ value somewhere between m^ and m 2 . The 
Incoherent Model estimates moistures very close to the 0~1 cm estimator. At the 
lower frequency m. values are slightly higher than m^ but as frequency is increased 
m. approaches m^ and at the highest frequency m. is slightly lower than m^ . The 
estimates from the Coherent Model are the lowest and predict erroneously zero 
moisture contents for a few profiles at 5.25 GHz and 7,25 GHz. There are also 
large variations in the soil moisture estimates. These variations are caused by the 
coherent nature of the model os was discussed earlier in section 3,3, 

Figure 10 shows soil moisture estimates os a function of frequency for profiles 
1, 9 and 2 which were considered dry, moist and very moist. For profiles 1 and 2 
the five estimates are quite similar while for the medium moist condition, the skin 
depth. Coherent and Incoherent Model give widely different estimates. 
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Depth, cm Depth, cm 


Soil Moisture, g/cm^ Soil Moisture, g/cm^ 



Soil Moisture, g/cm^ 

0.0 0.1 0.2 0.3 0.4 



Soil Moisture, g/cm^ 

0.0 0.1 0.2 0.3 0.4 



Figure 9. Soil moisture profiles for medium rough field. The 
profile number is indicated next to each profile. 
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Based on the above comparison of the various methods for representing a given 
moisture profile with o single moisture content value and based on the physical 
principles underlying each method, the Incoherent Model estimate m. has been chosen 
os the most appropriate. Hence, in all subsequent discussions, m. will be used 
exclusively. 

3.6 Spotiol Voriotions of Soil Moisture 

As was mentioned earlier in section 2, simultaneous with the radar measure” 
ments of a given field, soil samples were collected from each of eight samplir>g 
locations (Figure 1). After calculating m. fc<r each of the eight locations, the mean 
and standard deviation were obtained. Figure 11a**c show versus time graphs 
for the smooth, medium and rough fields for 0°, 10° and 20° incidence angles. Plotted 
on each figure are also soil moisture intervals (overage + standard deviation). The 
spatial variation of the soil moisture is appreciable and Is attributed mainly to the 
sprinkler irrigation used and to sampling uncertainty. The <^° response generally 
follows the soil moisture char>ges. Precipitation, which is shown on the graph, caused 
the surface roughness to change. For example in Figure 1 la, o° at nadir is larger 
for the data set taken on the 16th of July as compared to the 14th of July measure” 
ment while the soil moisture is actually lower. This is in consequence to the rain on 
the 14th which smoothed the field surface causing to increase. 


4.0 RESULTS 

A$ surface roughness plays an important role in the backscattering characteristics 
of bare surfaces, this section is devoted to discussiorw of the effect of roughness on 
the ongulor, spectral and moisture responses of o®. 

4.1 Angulor Response 

Figure 12 shows angular responses for the three bare fields with high levels of 
moisture content. 2.75 GHz, 5.25 GHz and 7.25 GHz were picked os representative 
frequencies for the low, mid and high range of the 2”8 GHz bond. The soil moisture 
content varies os a function of frequency for each of the surfoce roughnesses. This 
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Figure 12. Angular response of the scattering coefficient 

lor the smooth, medium rough, and rough fields 
for high levels cl moisture content at (al 2. 75 GHz. 
ibiS.KGHz, and (cl 7. 25 GHz. 
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does nof ind^ca^e different measuremenfs buf is inherent' in the moisture calculation 
of m.. As the frequency is increased the fields appear electromagnetically rougher. 

The combined angular response for all thiee surface roughnesses indicate that the 
effect of surface roughness is minimum (where the curves intersect) at approximately 
4® for 2.75 GHz, 10° for 5.25 GHz and 20° for 7.25 GHz. Do these curves also 
intersect at the same points under different moisture conditions? The answer is 
deferred to section 5. 

4.2 Spectral Response 

Figure 13 shows spectral responses for the three wetest field conditions for 0°, 
10° and 20°. At nadir (Figure 13a), between 5.25 GHz and 7.25 GHz cr° shows 
little variation, but between 2,75 GHz and 5.25 GHz fields M and R exhibit a 
decreasing response with frequency. <^° of field M, for example, decreases from 
9.9dBat 2.75 GHz to -1.5dB at 5.25 GHz. 

Figure 13b shows the spectral responses at 10°. a° of field R continues to 
decrease with frequency between 2,75 GHz and 5.25 GHz although at a slower rate, 
cf° of field M appears almost independent of frequency and field S indicares a strong 
increasing trend over the entire 2~8 GHz band. 

At 20°, unlike fields M and R which show only a weak response to frequency 
above 5.25 GHz (Figure 13c), ^° of field S increases rapidly. 

4.3 Soil Moisture Response 

The response of ct° to soil moisture at 2.75 GHz is shown in Figure 14a and 
14b for nadir and 10°, respectively. In each figure the data acquired from the three 
fields is indicated using different symbols. Linear regression lines are also shown 
along with the calculated correlation coefficient R and slope S of each. Since the 
slope is an indicator of the response of < 7 ° to moisture, it will be referred to as 
"sensitivity" S in dB/.Ol g/cm . The linear correlation coefficient and sensitivity 
calculated on the basis of all data points are also noted. 

At nadir (Figure 14a), the sensitivity of the smooth field is highest and the 
sensitivity of the rough field S^ is the smallest while at 10° (Figure 14b), the reverse 
is true. This reversal suggests that the effect of roughness may be minimized by operat- 
ing at an appropriate angle between 0° and 10°. The dramatic loss in sensitivity of the 
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Figure 13. Spectral response of the scattering coefficient for 
smooth, medium rough, and rough fields for high 
level of moisture content at lai 0° Inadir). Ibi 10®, 
and Ici 20° angle of incidence. 
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Figure 14. Soil moisture responses for the three surface roughness profiles at 2. 75 GHz for la) 0® (nadir) and (b) 10“ angle of incidence 
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smoofh field as fhe angle was increased from 0° fo 10° is affribufed to fhe smoofh 
appearance of fhe surface af 2.75 GHz. Since no radar data were acquired befween 
0° and 10°, an alfernafive fo invesfigaMng fhe moisfure response af an angle befween 
0° and 10° is fo invesfigafe fhe response af 10° buf af a higher frequency. Increasing 
fhe frequency causes fhe surfaces fo appear rougher v/hich is somewhaf comparable fo 
decreasing fhe angle of incidence. Figure 15 shows fhe moisfure response of 10° and 
4.75 GHz. The correlafion coefficienf calculafed for all fhe dafa poinfs regardless 
of surface roughness is 0.75. Considering fhe uncerfoinfy associafed wifh fhe measure 
menf of o° (Table 2) and fhe uncerfainfy associafed wifh fhe measuremenf of m. 
(Figure 11), fhe response shown in Figure 15 is very encouraging. Wifh a sensifivify 

of 0.25 dB/.Ol g/cm , fhe expecfed change in refurn befween moisfure levels of 

3 3 

0.05 g^cm (relafively dry) and 0.35 g/cm (very wef) is 7.5 dB . 


5.0 CHOICE OF OPTIMUM SYSTEM PARAMETERS 


In fhe previous secfion angular, specfral and soil moisfure responses were sfudied 
for all fhree roughnesses as a funcfion of fhe various sysfem paramefers. If was demon” 
sfrafed fhaf fhe effecfs of roughness can be reduced by proper choice of frequency and 
angle of incidence, while of fhe same fime refaining good sensifivify fosoil moisfure. 

In fhis secfion a more comprehensive analysis of fhe dependence of cr® on moisfure, 
surface roughness, frequency, polarizafion and angle of incidence is presenfed in a 
unified quanfifofive approach. 

5 . 1 Effecf of Roughness 

The paramefer chosen fo represenf roughness in fhis sfudy is rms heighf, h. 

As fhe radar responds fo bofh large scale and small scale roughness fhis represenfafion 
is undoubfedly nof opfimal. However as an inifial esfimafor (in ferms of ease of 
measuremenf) rms heighf seems fo be an adequafe quanfifier. Linear regression lines 
of cr° as a funcfion of soil moisfure confenf were calculafed for each individual surface 
roughness (rms heighfs of 0.88 cm, 2.6 and 4.3 cm) af each angle of incidence, 
polarizafion and frequency combinafion. The a° values shown in Figures 16-18 were 
calculafed using fhe regression equafions and selecfed values of soil moisfure. Figure 
16 shows o° versus h curves af nadir for 2.75 GHz, 4.75 GHz and 7.25 GHz. As 
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Frequency (GHz): 4.75 
Angle of Incidence: 10° 

Polarization: VV 

• Smooth Surface (RMS Height • 0.88 cm) 

▼ Medium Rough Surface (RMS Height = 2.6 cm) 
■ Rough Surface (RMS Height ■ 4.3 cm) 

S = Soil Moisture Sensitivity (dB/0.01 g/cm^) 

R ■ Correlation Coefficient 



Figure 15. Scattering coefficient response as a function of soil 
moisture for the combination of all three surface 
roughness profiles. 
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frequency Is Increased (a) the radar gets more sensitive to roughness for all soil 
moisture conditions and (b) sensitivity to soil moisture changes decreases. At 
2.75 GHz the radar Is Insensitive to surface roughness changes for values of h 
greater than 2.5 cm. Due to the large change of with h particularly for the 
0.35g/cm^ (14 dB at 2.75 GHz, 18 dB at 4.75 GHz and 16 dB at 7.25 GHz), 
operating at nadir would be an unsuitable choice. 

Figure 17 shows roughness responses ot four frequencies at 10°. As frequency 
Is Increased between 2.75 GHz (Figure 17a) and 7.25 GHz (Figure 17d), the total 
variation between a° of field S and a° of field" R decreases. Ideally, the optimum 
condition Is one where the moisture curves ore parallel horizontal lines. In 
addition to Independence to roughness, however. It Is desired that the separation 
between the <7° lines corresponding to different moisture levels be as far apart as 
possible. At 10°, the data shown In Figure 17c at 4.75 GHz, W polarization, 
satisfies these conditions best over the 2-8 GHz band. 

At 20° (Figure 18), variations In due to variations In h also decrease 
with frequency, but It appears that the condition of almost roughness Independence 
requires a frequency higher than 7.25 GHz. Increasing frequency or angle, 
however, results In a reduction In the magnitude of the sensitivity S (represented 
In the figures In terms of the separation between the different curves). 

5.2 Constont p° Contours 

Figures 16-18 have Indicated the o° response to rms height as a function 
of angle of Incidence and frequency. It Is also evident that for o given set of 
sensor parameters there exists more than one combination of moisture and rms height 
for which is fHe some. For an operational system, the soil moisture content Is 
to be predicted on the basis of a° measurements. To Investigate the degree of 
uncertainty associated with the prediction of mj from a°. Figures 19-21 were 
constructed. Using the curves shown In FIguresl6-18 , a family of constant 
contours were generoted with h and m. as axes for each frequencyangle combination. 

For minimum uncertainty in the estimate of m., which implies independence 
of h, the <^° contours should be as close to verticol lines as possible. Among the 
frequency-angle combinations shown in Figures 19-21, the data at 10° and 4.75 GHz 
(Figure 20c) approaches the desired condition the best. Since the data In Figures 
19-21 are the same data plotted earlier In Figures 16-18, the conclusion reached in 
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this secfion is fhe same os the conclusion reached in the previous section. The 
contour type of presentation is useful in determining the expected uncert iinty in 
the estimate of m. over o given ronge of h. For example, ol 4.75 GHz, 10°, and 
VV polarization, the m. ranges associated with the contours of Figure 20c over o 
range of h between 1 .0 cm and 4,0 cm ore os given in Table 3: 

Table 3: Mc^isture ranges corresponding to specific a° values at 4.75 GHz, 

10 ^ VV polarization with h varying between 1 .0 cm and 4,0 cm. 


o° (dB) 

m. (g/cm^) 

-10 

<.075 

- 7.5 

0.11 < m. < 0.17 

- 5.0 

0.24 

- 2.5 

0,3 < m. < 0.34 

0.0 

>0.36 


The contours depicted in Figures 19-21 suggest that the moisture estimate 
con perhaps be improved even further through the use of more than one frequency. 

Tills topic is, however, deferred to o future reporting, 

5.3 Regression Analysis 

The analyses presented in the previous sections indicate thot the angular 
region of interest in terms of soil moisture mapping Is between 0° (nadir) and 
20°. The radar data were acquired at five angles, 0° to 40° in 10° steps. In 
order to expand the number of available angles between 0° and 20°, non-linear 
interpolation of the angulor response was used to evaluate the values of o at 
3.3°, 6.7°, 13.3° and 16.7°. This procedure was repeoted for each frequency- 
polarization data set. Thus now . the new data base consists of eight frequencies 
X nine angles x two polorizotion = 144 data sets or <J° os a function of moisture 
for each of the three fields. Linear regression analysis wos performed for eoch of 
these data sets. After conducting this process for each of the moisture representations 
defirred in section 3, it was determined that the correlation coefficient was the 
highest for the overwhelming nrxjjority of the coses whet: moisture content was 
represented by m., the Incoherent Model definition. This result confirms its superior 
physical basis (in comporison to the other definitions of inoishi'.e content) in terms 
of the signal-target interaction process. 
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An example is shown in Figure 22 where the calculafed correlafion coefficient’ 

R between rf° and moisture content is plotted os o function of frequency for moisture 
content in the 0“1 cm layer, moisture content based on the Skin Depth model and moisture 
content based on the Incoherent Model. The results clearly indicate that overall, 

R is highest when mj is used to represent moisture content. 

The results of the regression analyses will be presented for three coses. First, 
results of the response to moisture for each individual field will be presented. 

Second, results of the cr° response to moisture for fields M and R together will be 
presented. This case was included to demonstrate the effect of adding field S, which 
is presented last. Since field S is almost an extreme case in terms of surface rough" 
ness of fields under natural conditions, the design of an operational radar for mapping 
so’l moisture should be based on an extrapolation between the results of the last two 
cases . 

For each of the three cases the procedure outlined below was followed: 
correlation coefficients and sensitivities were calculated for each angle of incidence, 
polarization and frequency combination. Then for a given angle of incidence and 
polarization combination the optimum (highest) correlation coefficient was picked. 

The frequency at which this optimum value occurred was noted along with the corre” 
spending sensitivity. To facilitate brevity the frequency and sensitivity associated 
with the optimum correlation coefficient are henceforth termed "optimum frequency" 
and "optimum sensitivity" respectively. 

5.3.1 Moisture Response for Individual Roughnesses 

Figure 23 presents the "optimum" radar perfonnance of field S as a function of 
angle of incidence. The optimum performance is defined in terms of three parameters: 
optimum frequency, optimum correlation coefficient and optimum sensitivity. It should 
be emphasized that at a given angle and for a given polarization, these three parameters 
are coupled together and hence, should be interpreted as a group. For example, the 
optimum correlation coefficients at 10° and 20° (Figure 23a) are approximately the 
same, and similarly are the optimum frequencies (Figure 23c), but the sensitivity at 
10° (Figure 23b) is larger than the sensitivity at 20° by about a factor of 1 .5. Hence, 
better performance would be expected at 10° compared with 20°. 
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Figure 22. Linear correlation coefficient between cr® and moisture content (with the 
latter expressed in terms of m; , m^, and mg) as a function of 
frequency. Data from all three fields were used in the regression analysis. 
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An Important note should be made regarding •'hese curves. Since the data used 
to generate the curves covers the band 2.75-7.25 GHz, the true optimum performance 
at some angles might be at frequencies outside this band. Consider, for example, the 
end points in Figure 23c. The shape of the curve suggests that at 0° the optimum frequency 
might have been lower than 2.75 GHz if such data were available. Similarly, at 40°, 
a frequency higher than 7.25 GHz might g've a higher correlation coefficient. In 
between these angles, howeve.', it is clear that a frequency around 5 GHz is optimum. 

Figures 24 and 25 present optimum performance curves for fields M and R respectively. 
Whereas the optimum frequency increases shnrpiv between 0° and 10° for field S (Figure 
23c), the optimum frequency curves for fiela R (Figure 25c) show a decreasing trend. 

5.3.2 Moisture Response for Medium and Ro ug h Fields Combined 

Optimum performance curves generuled on the basis of the combined data from 
fields M and R are shown in Figure 26. The optimum correlation is a maximum (0.89) 
at 10° and 2.75 GHz, for both polori? iti-ins. In the 0°-20° angular range, the optimum 
frequency varies within a narrow range between 2.75 GHz and 3.25 GHz. 

5.3.3 Moisture Response for Smooth, Medium and Rough Fields Combined 

When data from all three roughnesses are combined together, the optimum 
parameters show a strong dependence on angle of incidence. The optimum correlation 
coefficient (Figure 27a) increases from about 0.66 to 0.75 (for VV polarization) as 
the angle is varied between nadir and 10°, and then decreases rcpidly with angle. The 
corresponding sensitivity (Figure 27b) decreases approximately expe lentially with angle 
of incidence and the optimum frequency (Figure 27c) increases with angle from 3.25 GHz 
at nadir to 7.25 GHz for angles higher than 20°. The trend of the optimum frequency 
curves between nadir and 20° suggests that the optimum frequencies at angles higher 
than 20° would be higher than 7.25 GHz if such data were available . 
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6,0 RECOMMENDATIONS FOR AN OPERATIONAL SYSTEM 


A V ‘r earlier, the surface of field S Is a rather rare case In terms of 

tl ' maj irlty oi tef/t-*'' surfaces. Hence, by Including the data of field S to generate 
the perfoi xnce c ^f Figure 27, we are applying more stringent requirements 
on the per* • o* a ladcr s;.'stem for detecting soil moisture content than Is necessary 

■,r*der most '...l*’on« E. rluding field S, on the other hand, means a bias towords 
ro-.'/her surface FIguri Hence design parameters for an operational system can 

oe o fined as a a mpromise b<5*v/een those recommended by Figure 26 (Fields M and R 
V, iiy) and those recornmende'^ by Ftgu'® 27 (all fields). 

’‘.used on optimum correlotlon coe fficient and optimum sensitivity considerations, 
both figures Indicate that thr nguhr range extending between about 7° and 15° Is 
optimum. The corresponding optimu •' frequer. y, however. Is around 3 GHz In 
Figure 26 and around 5 GHz In FIgur. “^7, As a compromise, 4 GHz Is recommended. 
Over the 7°-15° < >aular range, ' ‘H ar \ VV p iJarIzatlons appear to have approximately 
the same soil mos/’T. i« "^ponse. 

In summaiy, the ri 'ommended roc..,r '-r.’'t-ieters are: 

Ar .> ',t Incidence Range: 7 -;J° 

Frequency; 4 GHz 
Polarizr . V HhorVV. 

It is most Interest' ig ‘.I j' rh above sensor 'rs recommended on the basis 

of radar measurements of gr,-uod agree perfectly wel* with 4-8 GHz measurements 
of vegetated fields reported by ^loby [5] . Data acquired for a variety of crop types were 
used to calculate the correlation coefficient between and soil molstu’e content. The 
highest correlation observed was at 10° and 4.7 GHz! 
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